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Abstract 

Residual  stresses  in  a  stack  of  the  anode- supported  planar  SOFC  were  measured  by  the  X-ray  diffraction  method  (the  sin2 method). 
The  stack  used  for  the  stress  measurements  was  composed  of  a  single  cell  and  separators  of  an  alloy.  In  this  measurement,  the  residual 
stresses  in  the  electrolyte  under  the  alloy  separator  and  the  cathode  were  focused  on.  In  order  to  detect  the  diffraction  from  the  electrolyte 
under  the  separator  or  the  cathode,  an  X-ray  energy  of  38.8  keV  was  selected.  For  the  stress  measurement,  a  diffraction  peak  of  YSZ(7  11) 
plane  was  used.  A  synchrotron  radiation  was  employed  as  an  excellent  X-ray  radiation  for  precise  stress  measurements.  In  addition  to  the 
stress  measurements,  numerical  simulations  for  the  residual  stresses  in  the  cell  stack  were  carried  out. 
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1.  Introduction 

Anode- supported  cells  are  inherently  strong  against 
stresses  because  the  electrolyte  hardly  suffers  a  tensile 
stress.  For  cell  stacks,  certain  alloys  can  be  used  as  inter¬ 
connectors  and  auxiliary  components  of  cell  stacks,  and 
thus,  using  the  flexibility  of  the  alloys,  it  is  possible  to  make 
the  cell  stack  robust  against  a  mechanical  load  [1,2].  On  the 
other  hand,  regarding  thermal  stresses,  great  care  should  be 
taken  into  account,  since  the  thermal  expansion  behavior  is 
different  between  the  cell  and  the  alloys.  The  thermal  stress 
induces  a  deformation  of  the  cell  stack,  and  sometimes,  the 
cell  stack  is  destroy.  To  reduce  the  thermal  stress  in  the  cell 
stack,  the  metal  alloys  have  to  be  selected  carefully  so  that 
the  thermal  expansion  coefficients  (TEC)  would  match  well 
with  that  of  the  cell.  To  study  the  best  alloy  as  the  stack 
component,  we  simulated  the  residual  stresses  in  the  cell 
stack  with  alloys  of  various  TEC.  However,  it  was  difficult 
to  calculate  the  residual  stress  in  the  cell  components  with  a 
high  accuracy  since  there  were  many  uncertain  parameters 
in  the  stress  calculation  (e.g.  a  temperature  at  which  the  cell 
components  are  constrained  each  other,  or  effect  of  the  re¬ 
duction  of  the  anode  on  the  stress).  To  evaluate  the  residual 
stress  in  the  cell  components  precisely  by  the  numerical  cal¬ 
culation,  it  is  indispensable  to  measure  the  residual  stresses 
practically,  and  fix  these  uncertain  parameters. 
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The  X-ray  stress  measurement  is  an  effective  method  to 
measure  the  residual  stress  in  materials  without  destruc¬ 
tion  [3-6].  We  applied  the  method  to  measure  the  resid¬ 
ual  stresses  in  the  cell  stack.  A  synchrotron  radiation  was 
used  as  an  excellent  X-ray  radiation  for  precise  stress  mea¬ 
surements.  Residual  stresses  in  the  electrolyte  of  the  cell 
stack  used  for  practical  measurements  of  the  stack  perfor¬ 
mance  were  measured  directly  using  the  X-ray.  In  addition 
to  the  naked  electrolyte,  the  stresses  in  the  electrolyte  under 
the  cathode  or  the  alloy  were  measured.  Numerical  calcu¬ 
lations  on  the  residual  stresses  in  the  cell  stack  were  also 
carried  out,  and  the  experimental  and  numerical  results  were 
compared. 


2.  Experimental 

Anode- supported  cells  were  fabricated  by  a  co-firing  pro¬ 
cess.  Details  of  the  cell  fabrication  have  been  reported  else¬ 
where  [7].  For  the  cell  stack,  electrically  conductive  or  insu¬ 
lating  alloys  were  applied  to  the  interconnectors  and  separa¬ 
tors,  respectively.  Fig.  1  shows  the  cross-sectional  diagram 
for  a  single-cell  stack.  A  cell  was  supported  using  a  thin  al¬ 
loy  foil,  and  the  alloy  foil  was  sandwiched  by  alloy  frames. 
A  glass  sealant  was  used  to  connect  the  cell  and  the  foil,  or 
the  foil  and  the  frame.  The  alloy  foil  has  a  role  to  separate 
the  fuel  and  air.  Metal  interconnectors  as  current  correctors 
were  set  at  the  upper  and  the  lower  side  of  the  cell,  but  they 
are  not  drawn  in  the  figure. 
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Fig.  1.  Cross-sectional  view  of  the  one-cell  stack. 


The  residual  stress  measurements  were  carried  out  using 
the  single-cell  stack.  The  cell  stack  has  already  been  tested 
its  performance  at  750  °C  in  an  electrical  furnace.  The  used 
stack  was  directly  mounted  on  the  sample  stage  at  the  syn¬ 
chrotron  facility  as  shown  in  Fig.  2  and  the  residual  stresses 
of  the  cell  stack  were  measured.  Fig.  3  shows  the  image 
photo  of  the  single-cell  stack.  The  center  black  square  is  a 
cathode,  and  around  the  cathode,  the  naked  electrolyte  ex¬ 
ists.  The  edge  part  of  the  electrolyte  is  covered  by  the  alloy 


Fig.  2.  Photo  image  of  the  experimental  stage  at  BL19B2  line. 


Fig.  3.  Photo  image  of  the  one-cell  stack.  The  incident  beam  was  irradiated  to  the  three  points  indicated  by  filled  circles. 
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foil,  and  the  foil  was  supported  by  the  separator  frame.  In 
this  stress  measurement,  we  focused  on  the  residual  stresses 
in  the  electrolyte  of  the  cell  stack  at  three  different  points. 
The  three  measuring  points  are  pointed  by  white  circles  in 
the  figure.  One  stress  focused  was  that  at  the  naked  elec¬ 
trolyte,  and  the  others  were  the  stresses  in  the  electrolytes 
under  the  cathode,  or  the  alloy  foil.  Although  the  residual 
stress  at  the  naked  electrolyte  can  be  measured  easily,  it  is 
difficult  to  measure  the  residual  stresses  in  the  electrolyte 
under  the  cathode,  or  the  alloy  foil.  To  measure  the  residual 
stress  in  the  electrolyte  under  the  cathode,  or  the  alloy  foil, 
it  is  necessary  to  detect  the  diffraction  of  the  X-ray  from  the 
electrolyte  through  the  cathode  or  the  alloy  foil.  Hence,  it 
was  important  to  select  an  appropriate  X-ray  energy  which 
enables  us  to  detect  the  diffraction  from  the  electrolyte  un¬ 
der  the  cathode  or  the  alloy.  Thicknesses  of  the  cathode  and 
the  foil  were  30  and  50  p,m,  respectively,  and,  the  shortest 
beam  path  in  the  alloy  was  100  jutm  even  when  the  beam  was 
irradiated  in  the  perpendicular  direction  to  the  cell  surface. 
The  main  element  of  the  alloy  was  a-Fe,  and  the  attenuation 
length  of  the  X-ray  in  the  iron  is  displayed  in  Fig.  4.  The 
X-ray  energy  of  more  than  30  keV  is  required  to  obtain  the 
attenuation  length  of  100  jam  in  a-Fe.  Thus,  we  selected  the 
X-ray  energy  of  38.7  keV  for  the  stress  measurement. 

A  synchrotron  radiation  was  employed  as  the  X-ray 
source.  The  synchrotron  radiation  is  an  ultra-bright,  highly 
directional,  and  collimated  light.  In  addition,  wide  spectrum 
of  wavelengths,  stretching  from  infrared  to  X-rays  can  be 
used.  For  the  present  stress  measurement,  an  X-ray  beam 
with  a  high  energy  and  a  high  brightness  are  necessary,  and 
thus,  the  synchrotron  radiation  is  the  best  source  for  the 
present  purpose. 

The  experiment  was  performed  at  BL19B2  line  at 
SPring- 8, which  is  a  medium-length  hard  X-ray  bending 
magnet  designed  for  the  purpose  of  engineering  science 
researches.  The  Huber  multi-axis  diffractometer  composed 


Table  1 

List  of  the  stress  measuring  conditions 

X-ray  source  Bending  magnet 

X-ray  energy  38.7  keV 

Monochrometer  Si3 1 1 

Measuring  method  Side-inclination 


of  conventional  four  circle  goniometer  and  additional  four 
axes  was  mounted  at  BL19B2,  which  was  suited  for  the 
present  stress  measurement. 

Table  1  shows  a  list  of  the  conditions  for  the  stress  mea¬ 
surements.  As  mentioned  above,  the  X-ray  energy  was  set  to 
38.7  keV,  which  was  enough  to  detect  the  diffraction  beam 
from  the  electrolyte  through  the  alloy  foil  or  the  cathode. 
Fig.  5  displays  the  result  of  0-20  scan  for  the  electrolyte  of 
the  cell  stack.  The  incident  beam  was  irradiated  to  the  naked 
electrolyte  part.  The  peaks  pointed  with  asterisks  correspond 
to  the  cubic  YSZ  structure,  and  the  other  peaks  are  diffrac¬ 
tions  from  Ni  of  the  Ni/YSZ  anode  under  the  electrolyte. 
For  the  X-ray  stress  measurement,  in  general,  a  diffraction 
peak  at  a  high  20  angle  is  used  because  enough  large  ^  can 
be  used.  On  the  other  hand,  the  intensities  of  the  diffraction 
peaks  at  a  high  20  angle  are  weaker,  and,  many  peaks  from 
different  diffracting  planes  overlap  each  other.  Thus  an  op¬ 
timal  diffracting  plane  has  to  be  selected  for  the  stress  mea¬ 
surement.  We  searched  a  single  diffraction  peak  existing  at 
a  higher  angle  and  having  enough  intensity  for  the  stress 
measurement,  and  selected  the  peak  of  YSZ(7  11)  plane.  As 
one  can  see  in  Fig.  5,  the  peak  center  of  YSZ(7  11)  peak 
is  around  25.7°.  If  the  iso-inclination  method  is  employed, 
available  maximum  \j/  becomes  <13°,  which  is  not  enough 
to  obtain  accurate  analyses.  To  obtain  an  enough  high  tilting 
angle  we  employed  side-inclination  method  instead  of  the 
iso-inclination  method.  Regarding  side-inclination  method, 
the  footprint  of  the  beam  at  the  surface  of  the  sample  at  a 


Fig.  4.  Attenuation  length  of  the  X-ray  inside  iron  when  the  direction  of  the  incident  beam  is  perpendicular  to  the  sample  surface. 
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Fig.  5.  X-ray  diffraction  pattern  for  the  electrolyte  of  the  anode-supported  cell. 


large  ifr  becomes  large.  The  size  of  the  incident  beam  at  the 
sample  surface  was  controlled  with  a  Ta-blade  slit  so  that 
the  incident  beam  would  be  irradiated  to  within  the  targeted 
area. 


3.  Results  and  discussion 

Fig.  6  shows  the  result  of  a  0-20  scan  where  the  incident 
beam  was  irradiated  to  the  foil  part.  One  can  see  that  the 
YSZ(7  11)  peak  appears  in  the  chart,  indicating  that  the 
diffraction  from  the  electrolyte  through  the  alloy  foil  could 
be  detected.  Fig.  7(a)-(c)  display  the  shift  of  the  YSZ(7  11) 


peak  with  t/r.  The  incident  beam  was  irradiated  to  the  naked 
electrolyte,  the  cathode  part,  and  the  alloy  foil  for  cases 
(a)-(c),  respectively.  For  all  the  cases,  the  YSZ(7  11)  peak 
shifts  toward  higher  angle  with  increasing  t/r.  Only  for  case 
(c),  can  one  see  the  split  of  the  YSZ(7  11)  peak.  We  consider 
that  the  split  may  come  from  an  imperfect  control  of  the 
incident  beam  size  (i.e.  the  beam  spot  may  slightly  spread  to 
the  naked  electrolyte  part).  The  peak  profiles  obtained  were 
fitted  using  the  Gaussian  function  as  shown  in  Fig.  8,  and 
the  peak  centers  were  determined.  Fig.  9(a)-(c)  show  the 
d- sin2t/r  diagrams  estimated  using  the  detected  YSZ(7  11) 
diffraction  peak.  Although  the  scattering  of  the  plot  is  larger 
for  case  (c),  d  decreases  linearly  with  sin2t/r  for  all  the  cases. 


20  (deg) 


Fig.  6.  Result  chart  of  6-26  scan.  The  incident  beam  was  irradiated  to  the  foil  part. 
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Fig.  7.  Results  of  26  scan  for  different  \[r.  The  incident  beam  was  irradiated 
to:  (a)  the  naked  electrolyte  part  of  the  foil  part;  (b)  the  cathode  part; 
and  (c)  the  alloy  foil  part. 

Thus,  we  can  estimate  the  residual  stresses  in  the  electrolyte 
using  the  following  formula. 

E  1  dd 

a  =  -  x  —  x  - - —  (1) 

1  +  v  do  dsinz\/f 

where  do  is  the  interplanar  spacing  under  a  stress-free  con¬ 
dition,  E  and  v  are  the  elastic  constant  and  Poisson’s  ratio. 
In  the  X-ray  stress  measurement,  E  and  v  must  be  measured 
using  the  X-ray  method.  In  the  present  stress  measurement, 
however,  we  did  not  measure  the  X-ray  elastic  constants. 


Fig.  8.  Fitting  for  the  obtained  peak  profiles  in  Fig.  7(c)  using  the  Gaussian 
function. 
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Fig.  9.  The  d-sin2^  diagram  measured  for  the  single  stack:  open  squares, 
open  circles,  and  open  triangles  are  for  the  naked  electrolyte,  the  elec¬ 
trolyte  under  the  cathode,  and  the  electrolyte  under  the  alloy  foil,  respec¬ 
tively. 

Thus  we  employed  the  calculated  value  of  the  elastic  con¬ 
stants  for  the  stress  estimation  (they  are  listed  in  Table  2). 
These  elastic  constants  were  calculated  from  the  reported 
elastic  constants  of  cubic  YSZ  [4,8]  using  the  Voight  model 
[9].  In  the  Voight  model,  the  elastic  constant  is  indepen¬ 
dent  of  the  diffracting  plane.  The  estimated  residual  stresses 
in  the  electrolyte  using  the  calculated  elastic  constants  are 


Table  2 

List  of  the  stiffness  of  the  cubic  YSZ  used  for  the  calculation  of  the 
elastic  constant 


Elastic  constant 

Electrolyte  (GPa) 

Cn 

400 

C12 

110 

C44 

55 

E/(  1  +  v) 

182 
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Table  3 


List  of  the  estimated  residual  stresses  in  the  electrolyte 


Part 

Evaluated  stress  (MPa) 

At  naked  electrolyte 

486 

Under  cathode 

471 

Under  alloy 

551 

listed  in  Table  3.  Everywhere  in  the  electrolyte,  the  estimated 
residual  stresses  were  compressive  stresses.  As  compared 
to  in  the  naked  electrolyte,  the  calculated  residual  stress  is 
slightly  small  under  the  cathode,  and  large  under  the  alloy 
foil.  Thus,  it  was  found  that  stacking  the  cell  changes  the 
residual  stress  in  the  electrolyte. 

In  addition  to  the  stress  measurements,  the  distributions 
of  the  residual  stresses  in  the  stack  were  simulated  us¬ 
ing  a  computer  simulation  technique.  The  stack  geometry 
(cathode-electrolyte-anode,  separator,  and  frame)  was  mod¬ 
eled  precisely,  and  the  distributions  of  the  residual  stresses 
in  the  stack  components  were  calculated  using  the  finite  ele¬ 
ment  method  (details  of  the  method  for  the  stress  simulation 
have  been  reported  previously  [10]).  In  the  present  simu¬ 
lation,  we  assumed  several  calculating  conditions.  First  as¬ 
sumption  was  the  temperature  at  which  the  stack  components 
were  constrained  (connecting  temperature).  The  tempera¬ 
tures  at  which  the  electrolyte  and  the  anode,  the  electrolyte 
and  the  cathode,  and  the  electrolyte  and  the  alloy  foil  are 
connected  were  assumed  to  be  1400,  1250,  and  900  °C,  re¬ 
spectively.  In  the  simulation,  above  1400,  1250,  and  900  °C, 
the  anode,  the  cathode,  and  the  glass  sealant  were  treated  as 
plastic  materials.  That  is,  only  below  the  connecting  tem¬ 
perature,  these  components  would  contribute  to  the  residual 
stress  in  the  electrolyte.  Second  assumption  was  that  the  re¬ 
duction  of  the  anode  (i.e.  the  change  of  NiO  to  Ni  in  the 
anode)  would  not  affect  the  residual  stress  in  the  electrolyte. 
Our  previous  stress  measurements  for  the  co-fired  cells  re¬ 
vealed  that  there  was  no  marked  difference  in  the  residual 
stresses  between  the  cells  with  the  oxidized  and  the  reduced 
anodes  [10]. 

Fig.  10(a)  displays  the  simulated  stress  distribution  in  the 
electrolyte  at  room  temperature.  The  absolute  value  of  the 
residual  stress  is  smaller  under  the  cathode  and  larger  un¬ 
der  the  alloy  foil  than  in  the  naked  electrolyte,  which  is 
qualitatively  consistent  with  the  result  of  the  stress  measure¬ 
ments.  The  comparison  of  the  residual  stresses  for  the  mea¬ 
sured  and  simulate  values  are  listed  in  Table  4.  As  compared 
to  the  simulated  values,  the  measured  values  are  larger  by 


Table  4 


Comparison  of  the  measured  and  simulated  stresses 


Part 

Measured 
stress  (MPa) 

Simulated 
stress  (MPa) 

Modified 

measured 
stress  (MPa) 

At  naked  electrolyte 

486 

679 

633 

Under  cathode 

471 

670 

612 

Under  alloy 

551 

705 

719 

Principal  Stress 
[MPa] 


max  = 


mill  =- 


>-6.60e+0: 

<-6.60e+0: 

<-6.66e+0i 

<-6.72e+0: 

<-6.78e+0: 

<-6.84e+0: 

<-6.90e+0; 

<-6.96e+0J 

-5.75e+02 

-6.96e+02 


(b) 


Fig.  10.  Simulated  distribution  of  the  residual  stress  in  the  electrolyte  at 
room  temperature:  (a)  for  the  single  stack;  and  (b)  for  the  single  cell. 


20-30%.  This  discrepancy  may  originate  from  the  follow¬ 
ing  two  reasons.  One  is  the  connecting  temperatures  were 
not  proper  in  the  stress  simulation.  If  the  actual  connecting 
temperature  is  lower  than  we  assumed  it  to  be,  the  simu¬ 
lated  values  become  larger  than  the  measured  values.  An¬ 
other  reason  is  that  the  elastic  constants  used  for  the  stress 
estimation  were  inadequate.  Since  the  X-ray  elastic  con¬ 
stants,  in  general,  depend  on  the  diffracting  plane,  also  in 
the  present  stress  measurement,  the  X-ray  elastic  constants 
must  be  measured  using  the  YSZ(7  11)  peak.  However,  we 
used  the  calculated  values  as  the  elastic  constant,  and  thus  it 
may  result  in  a  less  estimation  of  the  residual  stresses.  Fig. 
10(a)  and  (b)  are  the  comparison  of  the  simulated  residual 
stresses  in  the  electrolytes  for  the  single  cell  and  the  single 
stack.  From  the  figures,  we  can  find  that  the  residual  stress  at 
the  naked  electrolyte  for  the  stack  model  is  almost  the  same 
as  that  at  the  center  for  the  single  cell.  This  suggests  that  the 
stacking  of  the  cell  will  not  affect  the  residual  stress  at  the 
naked  electrolyte.  In  our  previous  evaluation  of  the  residual 
stress  for  the  single  cell,  the  measured  stress  and  the  simu¬ 
lated  stress  in  the  electrolyte  were  almost  the  same.  In  the 
previous  evaluation,  the  same  X-ray  elastic  constants  as  the 
present  estimation  were  employed  to  estimate  the  residual 
stress.  Besides,  in  the  previous  stress  calculation,  the  con¬ 
necting  temperature  between  the  electrolyte  and  the  anode 
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was  assumed  to  be  1400  °C  which  is  the  same  as  the  present 
assumed  values.  Considering  these  facts,  we  can  conclude 
that  the  origin  of  the  discrepancy  between  the  measured  and 
the  calculated  stresses  in  the  present  case  is  improper  elastic 
constants. 

As  mentioned  above,  for  the  stress  simulation,  the  resid¬ 
ual  stress  at  the  naked  electrolyte  of  the  stack  is  almost  the 
same  as  that  at  the  center  of  the  single  cell.  Accordingly,  we 
assume  that  the  residual  stress  in  the  naked  electrolyte  of  the 
stack  is  the  same  as  that  at  the  center  part  of  the  electrolyte 
for  the  single  cell,  and  modify  the  elastic  constants.  Table  4 
shows  the  modified  residual  stresses  in  the  electrolyte  as¬ 
suming  that  the  stress  at  the  naked  electrolyte  is  the  same  as 
that  of  the  single  cell.  The  modified  and  simulated  residual 
stresses  become  close.  For  more  accurate  experimental  esti¬ 
mations,  the  X-ray  elastic  constant  must  be  measured  using 
the  YSZ(7  1  1)  peak. 

The  strength  of  YSZ  against  the  compressive  stress  is 
more  than  1  GPa  [ll]and  thus  these  residual  stresses  are  not 
so  large  as  to  induce  a  failure  of  the  electrolyte. 

4.  Conclusion 

The  residual  stresses  in  the  electrolyte  of  the  single-cell 
stack  were  evaluated  using  experimental  and  numerical 
methods.  The  estimated  stresses  in  the  electrolyte  were 
compressive  stresses  of  633,  612  and  719  MPa  at  the  naked 
part,  under  the  cathode,  and  under  the  alloy  foil,  respec¬ 
tively.  Although  the  stacking  of  the  cell  slightly  increases 
the  residual  stress  in  the  electrolyte  under  the  alloy  foil,  the 


residual  stresses  are  not  so  large  as  to  induce  a  failure  of 
the  electrolyte. 
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